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Abstract 

The two-dimensional one-component logarithmic Coulomb gas is mapped onto a non- 
hermitian fermionic field theory. At (3 = 2, the field theory is free. Correlation functions 
are calculated and a perturbation theory is discussed for extending to other (3. A phase 
transition is found at the mean-field level at large f3. Some results are extended to spaces of 
constant negative curvature. 

1 Introduction 

The problem of the statistical mechanics of a system of particles in a two-dimensional plane, 
all possessing the same sign of charge and interacting via a repulsive Coulomb interaction, has 
been looked at by a number of authors]!]]. This is referred to as a one-component plasma since 
all particles have the same sign of charge. The Coulomb interaction is logarithmic in two- 
dimensions which makes it possible to introduce a dimensionless inverse temperature (5 such 
that the partition function of an TV-particle system is equal to 

N N N 

(n dzidzi) n \zi - Zj f n ^ 

i=l i<j i=l 

where U (r) is some background potential and where z = x + it and z = x — it are coordinates 
in the complex plane. 

It is believed that at some (3 of approximately 140 there is a first order phase transition 
in the system M. Below the phase transition the system acquires long-range orientational order 
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but does not acquire long-range positional order. The system is thus of interest as a model of 
crystallization. In the last section of the paper, we consider this model on a space of constant 
negative curvature. This may make it possible to look at effects of frustration since it is not 
possible to form a hexagonal crystal lattice without defects on a negative curvature space||. 

This problem is exactly solvable at j3 = 2 for a variety of background potentials. We 
introduce a new method of solving the problem, based on mapping to a non-hermitian field 
theory, which may be easier to deal with for certain potentials. In addition, we look at a 
perturbation theory which permits a mean-field analysis of the phase transition in this theory. 



2 Fermion Mapping at (3 = 2 

As a a statistical mechanics system, the model, for a finite number of particles in a uniform 
background, is defined by the partition function 

. N N N 

j (H d Zi dzt) n n - z 5 f n e -^ 2 (2) 

i=l i<j i=l 

N 

We note that the factor JT \z{— ZjY in the partition function is equivalent to a correlation 

i<j 

function in a free bosonic field theory: 

N 
i=l 

where the field $ has the action S = -a j dxdt(V<&) 2 . 

We may then go to an infinite number of particles, and write the partition function as 

Z = J [d®]e~ s+ 1 'e*^*^''' -i P V^${x,t) dx dt ^ 

It may be seen that, by perturbatively expanding the partition function in powers of e Jv/ ^"*, we 
recover the original statistical mechanics problem. 

This equivalence between eqs. (||) and (||) is similar to the equivalence between the 
statistical mechanics of an plasma consisting of both plus and minus charges and the field 
theory of the sine-Gordon equation. 

As is well known, the sine-Gordon equation at (3 = 2 maps onto a problem of free massive 
fermions, while at other temperatures the equation maps onto the Thirring model, a model of 
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interacting massive fermions We will follow an analogous procedure in this case, leading to 
a non-hermitian field theory of fermions. All correlation functions will be given at (3 = 2, while 
a perturbation theory in a four-fermion interaction will lead to other temperatures. 

Each term of the bosonic action translates into a given term of an action for relativistic 
fermions. The action S for the field $ translates into 2 J dx dttp R d^ipR + ijr L d z ij)L- The term 
e iV4n$(x,t) t rans i a tes into 27ra^ L il>ji, where a is an ultraviolet cutoff for the theory. The existence 
of only one kind of charge in the statistical mechanics theory leads to a non-hermitian field theory. 

It is seen that placing a charge in the statistical mechanics theory at a given point 
corresponds to turning a fermion at that point from a right-mover into a left-mover. The effect 
of the neutralizing background charge is to turn left-movers back into right-movers via the 
anomaly. One may imagine the neutralizing background charge as a charge at infinity. We can 
write i<3? as i&dft and integrate by parts to turn this into itdt& plus a boundary term. This 
then bosonizes into ^/7rtJ x , which is ^/Trt^^R — i/j^i/jr). This has the effect of introducing a 
time-dependent chemical potential which has different signs for the right- and left-movers. If 
we imagine following the Hamiltonian evolution of the fermionic field theory in imaginary time, 
right-moving particles are constantly converted into left moving particles, but the energy of the 
states keeps changing so that on average the number of holes in the negative energy sea is the 
same for the both chiralities. The effect of the boundary term resulting from the integration by 
parts is that we must start at t = — oo with a state in which all right-moving states (of both 
positive and negative energy) are filled and all left-moving states are empty and end at t = oo 
with a state in which all left-moving states are filled and all right-moving states are empty. Let 
us refer to the state at t = — oo as the state V~ . We will refer to the state at t = +00 as the 
state V + . 

The final action for the fermion field is 

J dx dt 2{i> R d^ R + ^\d t ij} L ) + V2^^y R + P t2^y L - rP R ip R ) (5) 

There are several ways of writing the background charge which correspond to choosing different 
gauges. The factor \f2p is chosen to cancel a factor arising later from an integral of J dte~ 2npt2 . 
It also serves to give the term in the action the correct dimensions. This factor is unimportant 
since the fixed density of the background charge means, due to charge neutrality, that the number 
of particles in the statistical mechanics problem is fixed and thus this factor just multiplies the 
partition function but does not affect the physics. 
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3 Correlation Functions at (3 = 2 



We compute correlation functions for the fermionic field and use them to then obtain correlation 
functions in the statistical mechanics problem. It will hopefully be clear in which theory a 
correlation function is being calculated. 

We will compute propagators for the fermion field as follows: we will first compute the 
propagators for the case in which the two operators lie on the line t = 0. Then we will, in 
an appropriate gauge for the background charge, generalize to arbitrary position of the two 
operators. 

In the operator formalism, we can write a two-point correlation function of fields iib{\) 1 ip{2) 
as the expectation value 

oo 

-i - J H{t)dt - J H(t)dt 

-{V + \e o ^(l)^(2)e — \V~) (6) 

Zj 



where 



- / H(t)dt - J H(t)dt 

Z = (V + \e o e -oo \ V -) (7) 



The term in the Hamiltonian, J dx dttpt tpR, can be written in Fourier components as 
J T^a) (k) La(—k) ft- Each term in the integral is an eigenoperator of the rest of the Hamiltonian, 
with eigenvalue 2k + Airpt, and so it is easy to rewrite eq. (o) in terms of operators acting just 
at time t = as follows 



(Je-^ +aW *)>/W(*)i«(-*)B / e- 2 * pt2+2kt dt)y%>aHk) L a{-k) R 

—{V + \e o V(l)^(2)e -°° \V~ ) 

Zj 

(8) 

The simplest two-point function to compute is (a^(A;)x,a(— k)n). In the expectation value 
of eq. (||), for each k the operators a){k)L and a(k)n must appear exactly once. It is clear then 

oo 

/ y^e~^P t2 + 2kt dt 

that Z = n e-°° . This is 

k 

fc 2 

Y[e^p (9) 

k 

By inserting the operator a' (fc)ia(— U)r into the expectation value, we remove one term from 

fc 2 

the product of eq. @. This term is e~p . Therefore, {a) (k) La{—k) r) = e ~p . In real space 
E e ~2P x ; and in an appropriate gauge it becomes \ ^&~ 2 ■ 
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We next compute the function (o^(A;)lo(— k)L- We note that initially all left-moving 
states are unoccupied at t = — oo. Then all states are occupied at t = oo. If the state with 
momentum k is occupied at t = then this gives a contribution of 1 to the desired two-point 
function. If it is unoccupied, then there is a contribution of 0. Therefore the desired two-point 
function is 



j e -2npt 2 +2kt dt 



J e -2npt^+2kt dt 



(10) 



In real space, this is n^T-re 71 2 W ; where an appropriate gauge must be used to extend correlation 



2niz 

functions off the line t = 0. 



The function in which both left-moving operators are replaced by right-moving operators 
is similar. The result is ^= e _7r ^' z ' . 

The most interesting two-point function to compute is (a"l"(fc)^a(— This requires 
that one a)(k)La{— k)ji operator appear at time t < and one must appear at time t > 0. 
Therefore, the two-point function is 

J e ~2Trpt 2 +2kt dt , J e -2npt'*+2kt> df 

V^-P- 56 " (11) 

J e -2«pf>+2kt dt 

— oo 

In real space, this is 2p j e^P^dt' Je-^dt J M e -£- pe k(ix+2t+2t>)_ That becomes 

-00 

00 

p j e-^dt 1 J e- 2 ^ t2 dte^ x+2t+2t '^^Tp (12) 



00 00 

This is equal to V2pl / dt' J dte^^ x2+m ' +4ix ^ t+t ">\ Integrate t' to obtain ±J\ J dt-^e^^+^l 
-00 n v + 2 



This is equal to ^\j\& r % x2 J^, x 2 e t /t which becomes 

i-y| e ^N 2 Ei ( vrp|z| 2 ) (13) 

where Ei is the exponential integral function. 

Having computed the propagators, it is trivial to compute correlation functions in the 
statistical mechanics theory. Let us compute the two-point particle correlation function at (3 = 2. 
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This is obtained by the expectation value (y/2~p^j^(0)^ R ,{0)\/^py[ j {z)'ipii(z)). There are two ways 
to pair off the operators in this expression using Wick's theorem. The result is 

p 2 (l - e-^M 2 ) (14) 

which agrees with the known result. 

4 Perturbation Theory 

By adding to the action of the fermionic theory a term of the form (tpTp) 2 , we will, due to the 
bosonization rules, change the value of (3. An attractive interaction will raise j3, while a repulsive 
interaction will lower (3. We will employ a simple mean-field type procedure: the interaction 
will be decoupled into an interaction with an external gauge field. Then, we will integrate out 
the fermionic field to obtain an action for the gauge field. It will be found that at sufficiently 
high, finite (3 the quadratic term in the action for the gauge field will change sign. This will be 
identified as the location of a second order phase transition. It will then be argued that cubic 
terms appear in the action for the gauge field, and will convert the phase transition to first 
order. 

For a given value of 0, the desired fermionic action is 

/l 1 
dx dt 2{p R d&i) R + r L 8 z ip L ) + v^p^y R + pt2ir{i>y L - ipy R ) + (- - - ^ ■. ^y R ^y L ■. 

(15) 

We can write the interaction term by introducing an external gauge field A R l- The action for 
the gauge field plus the coupling between the gauge field and the fermion field is 

/ dx dt A R J R + A L J L + -^-^-^ArAl (16) 

where we define J R =: ip R ip R : and Jl ='■ ipiipL '■ and where A R is the complex conjugate of Al- 
It will turn out, when we compute various diagrams, that results will not be gauge 
invariant. This is unrelated with the particle choice of gauge we made to compute propagators, 
and instead is a result of the boundary conditions due to the neutralizing background field. 
The theory is also not invariant under charge conjugation. This implies that Furry's theorem^] 
need not hold, and we may obtain contributions from fermion diagrams with an odd number of 
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photon vertices attached. This will then convert the second order transition into a weakly first 
order transition. 

Let us first compute the quadratic terms in the action for A R; l resulting from integrating 
out V- There are three different types of diagrams which must be considered, which involve 
computing different current-current correlation functions. We may have to compute a correlation 
function involving two currents of the same chirality. This will be considered first. We may have 
to compute a correlation function involving currents of the opposite chirality. This correlation 
function splits into two pieces: one piece which may be obtained by naively calculating diagrams, 
and one piece which is a result of the need for a regulator field to get rid of divergences. 

The introduction of the quadratic term will change the action of the gauge field, as a 
function of wavevector k, to 

3^TTxAr(&)Al(— - \{J R {k)J R (-k))A R (k)A R {-k) 

2 (17) 

-\{J L {k)J L (-k))A L {k)A L {-k) - (J R (k)J L (-k))A R (k)A L (-k)) 

All that needs to be computed now are some current-current correlators. 

The current-current correlation function {J R (0)J R (z)) is equal to ^nzp e " np ^ ■ We 
must take the Fourier transform of this at a given momentum k. In order to do this, it is 
necessary to introduce a massive regulator field. Without the term e~ wp \ z \ 2 in the correlator, 

9 i k 2 

the Fourier transform of 1/z is known to be j^jj^i, where k^ is equal to k x — ik t , kji is equal 
to k x + ikt, and \k\ 2 = kik R . The term e" 71 "^ 2 ' 2 multiplies the rest of the correlator in real 

space and thus convolves with the correlator in Fourier space. The k-th Fourier component of 

r i , l fc -'l 2 

the current-current correlator in momentum space is then J dk^dk R -^^e **p . This integral 

i k 2 Ifel 2 

can be performed analytically and the result is ^^(l — e ^p). It is seen that at large k 
the current-current correlator is unchanged from the correlator for a fermionic system with no 
non-hermitian term. The correlation function of two left-moving currents is calculated in the 
same way, simply replacing k^ by k R and vice- versa. 

The naive contribution to the correlator of two currents of opposite chirality is given 

oo 

by (J R (0)J L (z)). This is ^Ei(vrp|z| 2 ). The Fourier transform is ^Jdxdt J ^ e ^x-n pa (x +t )_ 

i 00 j |fc| 2 

Doing the integral over x and t first, we obtain J 4a7T p . Changing from a to 1/a, this 
equals ^(1 — e 4 *p). 

The necessity of introducing a massive regulator field to compute the correlation function 
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of two currents of the same chirality gives an additional contribution to the correlation function 
of two currents of opposite chirality. This piece is infinitely short range, and so independent 
of k. It is equal to For the fermionic system with no non-hermitian term, such a piece is 
needed to maintain gauge invariance. 

Even without looking closely at the Fourier transform of the current-current correlators, 
we may easily see that there will be an instability for finite (3 at the level of the quadratic action 
for Al^r. For the fermionic system with no non-hermitian term, the quadratic term changes 
sign, for all k transverse to the direction of the gauge field, at (3 = oo, which corresponds to a 
finite value of the attractive interaction. By adding the non-hermitian terms to the action, we 
have altered the quadratic terms. The correlation function of two currents of the same chirality 
has been reduced, but the reduction decays exponentially for large \k\ 2 . The correlation function 
of two currents of the opposite chirality has been increased by an amount which decays only 
algebraically for large \ k\ 2 . Thus, for sufficiently large \ k\ 2 , the term which is added to the action 
of the Al^r field is increased and the intability will occur at a lower value of the attractive 
interaction. This corresponds to an instability at a finite value of (3. 

Let us locate the transition temperature in this lowest order theory. The response to a 
gauge field is strongest when the field is transverse. The response then is (summing all different 
terms and taking Ar = Al) 

1 1 __fci 1 

— + — (1-e 4*,)+ (i_ e *" P ) (18) 

47T 47T k 2 

This function is equal to ^ at k = and k = oo. The function increases in magnitude as k 
increases from 0, until it hits a maximum, then decreases again. Numerically, the maximum is 
at k = 4.74711^. At this wavevector, we find that the theory goes unstable at a temperature 
of = 15.4036. 

In a mean-field approximation, in the absence of cubic terms, the instability discussed 
in the quadratic action for Ar^, would lead to a second order transition. It will now be shown 
that cubic terms are non-vanishing and that this transition becomes first order. Due to the 
complexity of the diagrams, we will not explicitly compute the cubic and quartic terms, and will 
simply show that the cubic term is non-vanishing. 

In 2 dimensions Furry's theorem is very easy to prove for a massive or massless Dirac 
field. If we have a fermion loop with an odd number of vertices, we may imagine reversing the 
direction of the propagation of the fermion around the loop. This will change the sign of all 
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propagators which preserve chirality (propagators of the form (i/Jri/jr) or (ip L ipL)) of the fermion 
field, and preserve the sign of all propagators that change chirality. Since the fermion must have 
the same chirality after it goes around the loop, there are an even number of propagators which 
change chirality, and thus an odd number of propagators which preserve chirality. This means 
that the total sign of the diagram changes when reversing the direction of propagation around 
the loop and thus the total contribution from both directions is zero. 

For the non-hermitian theory considered here, the above proof breaks down. When 
reversing the direction of propagation around the loop, we may change some propagators from 
(ip\ i {zi)tpR(z2)) to (ip\ i (z2)i^ 'l(zi)} ■ Since the non-hermitian theory lacks charge conjugation 
invariance, these two propagators have different values, and thus the two contributions do not 
cancel. This means that cubic terms do appear in the action for Ar^ and the transition becomes 
first order. It should be noted that these cubic terms are intrinsically not gauge invariant, since 
they vanish in a gauge in which Al = or Ar = 0. 



5 Extension to Curved Space 

We will also consider this model on a space of constant negative curvature. We may define such 
a space in the half-plane given by Im(z) > 0, or equivalently, t > 0. We will use a conformal 
gauge for the metric such that ds 2 = ji{dx 2 + dt 2 ). Then the scalar curvature R = — 1. 

It is easy to transcribe the action for the bosonic field to this curved space and the result 

is 

Z = J [d^S+jzJe^^-^pV^^dxdt (19) 

When going to a fermionic action, one must be slightly careful. The curved space means 
that in order to maintain the same ultraviolet cutoff a when measured in length Vds 2 every- 
where, the ultraviolet cutoff when measured in length \J dx 2 + dt 2 must vary proportional to t. 
Therefore, the operator J e l ^ r ^^{ x > t ) bosonizes into something proportional to {2irt)il^ L '^R 

This means that the final fermionic action for the curved space problem is 

J dxdt2(i/j R d^ R + + + (V&^fl " A^>l) (20) 

In this action, the density, p, is dimensionless, as we have picked units in which the curvature 
is one. We can use 1/p to measure the curvature of the manifold in units of particle spacing 
squared. As p — > oo, we expect to recover the flat space results. 
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It is now possible to proceed as before and calculate propagators. We will simplify and 
sketch the calculation of only (tp^tpji). We restrict to the line t = 1. The only new feature that 
emerges is that the final state is different. The initial state, |^ _ ), is still a state in which all 
right-moving states are occupied and all left-moving states are empty. The system starts in this 
state at time t = 0. However, there is no term left from the integration by parts at t = oo and 
therefore the system ends in the vacuum state \ V) which is the normal vacuum state for a the 
free fermion theory, where all negative energy states are filled and all positive energy states are 
empty. 

Following a similar procedure as before, we move all i/j^i/jr operators resulting from the 
Hamiltonian evolution to the time t = 1. There, an operator (k)La(—k)R comes with the 
amplitude 

i 

°? l - J 2k+8np/t'dt' 

Jdt-e « (21) 
o 

oo 

This is equal to / dt ^X^p e 2fc ( t ~ 1 ) , which is proportional to e~ 2k , , . Therefore, (at (k) ia{—k) r) 

o 1 1 

at t = is proportional to e 2k \k\ SlTp . Due to the initial and final states, \V~) and \V), the prop- 
agator is non- vanishing only for k < 0, when this Fourier transforms to .J; 1+g . 

The result for the correlation function in the statistical mechanics problem, obtained 
from calculating (0)liP(0)riP(x)liP(x)r), is 

1 



(rr^) 1+87Tp (22) 

4 



If we take a limit as p goes to infinity, and appropriately rescale x to measure length in the 
local particle spacing, the correlation function turns back into a Gaussian, as in the flat space 
case. Recall that the distance between two points on the line t = 1, located at x = and 
x = x is, measured in the metric for the negative curvature space, less than x. In reality it goes 
logarithmically with x for large x. Thus, the actual decay of correlations in the curved space 
problem is exponential, as would be expected from a high-temperature expoansion. The slower 
decay of the correlations in the curved space case, exponential instead of Gaussian, indicates 
that screening is less effective than in flat space. 
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6 Conclusion 

A new formulation has been given of the one-component plasma problem. This provides an 
alternative way of deriving old results. This technique may turn out to be simpler than others 
for certain background potentials. In addition, the perturbation theory of this model in terms 
of fermionic operators is very different from the previously developed perturbation theory for 
the statistical mechanics theory ||, and may lead to easier calculations. 

In particular, we located an instability of the theory at the quadratic level. It must 
be noted that the location of the second order transition in this theory, at j3 = 15.4036, is 
approximately an order of magnitude lower than the actual location of a first order transition. 
Also, the appearance of cubic terms, while converting the second order transition to first order, 
will further lower the transition (3. This is evidence that the higher order corrections to the 
effective action for the gauge field are not negligible. However, the manner in which the transition 
temperature is calculated makes this number very sensitive to small adjustments in the effective 
action; one must calculate the response functions in the fermionic theory and then relate (3 to the 
reciprocal of the change in the action for the gauge theory. The act of taking the reciprocal makes 
this procedure less accurate. Unfortunately, we are also unable to understand the particular 
wavevector which goes unstable, as we do not see how to easily relate this wavevector to any 
spacing in a triangular lattice. In addition, since the order parameter is a vector, it is possible to 
construct a mean-field state that breaks orientational symmetry but not translational symmetry 
(constant non- vanishing gauge field). This would require that the first wavevector to go unstable 
would be at k = 0. This does not happen yet at the one-loop level, though one would expect it 
would happen to higher orders. 

However, we may hope that the perturbation series for the effective action of the gauge 
field will be convergent, order by order in (1/2 — 1/(3), as the only actual instability of the theory 
is at (3 = oo. This implies a radius of convergence of 1/2, and thus the theory should converge 
for 1 < (3 < oo. Arguments like Dyson's instability argument for the non-convergence of the 
perturbation series do not apply here since there is no instability near p = 2@. To make this 
convergence more clear, it may help to adjust units so that the action for the gauge field is 



to avoid what may look like a singularity at (3 = 2. Thus, a sufficiently high order calculation in 




(23) 
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this theory should yield the effective action for the gauge field, which can then be treated in a 
mean-field fashion to, in principle, extract the transition properties with arbitrary accuracy. Of 
course, the series might diverge at the point of a second order phase transition due to infrared 
problems, but since a first order phase transition is expected to occur before the second order 
transition, this is not a problem. 

Finally, the model has been considered on a curved manifold. The two-point correlations 
in the original statistical model have an exponential decay instead of a Gaussian decay. Since 
the decay of this quantity reflects the effects of screening, it seems that on a curved manifold the 
system screens less well. This may be due to an effect of frustration, introduced by the curvature. 
It would be interesting to extend the perturbation theory to a curved manifold, both to look 
at the RPA as well as to look at correlation functions away from (5 = 2. It is believed from a 
perturbation theory for the original statistical mechanics model H that the correlation functions 
on a flat space begin to show short-range order as soon as (3 > 2. This means a lowest order 
perturbation theory calculation for the curved space may show interesting frustration effects. 
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